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Abstract 

We analyze an initial-boundary value problem for the Ostrovsky- 
Vakhnenko equation 
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on the half-line. This equation can be viewed as the short wave model 
for the Degasperis-Procesi (DP) equation. We show that the solution 
u{x, t) can be recovered from its initial and boundary values via the 
solution of a 3 X 3 vector Riemann-Hilbert problem formulated in the 
complex plane of a spectral parameter z. 
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1 Introduction 


In this paper, we consider the partial differential equation (PDEs) 

Uxxt — “ibUx + ‘iUxUxx + UUxxx = 0 (1-1) 

where 6 > 0 is a parameter and u = u{x, t) is real-valued. This equation 
stems from the short-wave limit of the Degasperis-Procesi equation (DP) |T], 
which is a model of nonlinear shallow water waves: 

Ut Uxxt T ^(XUx 4:UUx - ^UxUxX T ^^XXX- (1'^) 

Indeed, introducing new space-time variables and a scaling of u by 

, X ! I 

X — —, t — U = —;r, 

where £ is a small positive parameter, then fll.ip is the leading term of fll.2p 
as £ —)■ 0. Thus, equation fll.ll) can be named as the short wave model of the 
Degasperis-Procesi equation. 

Equation fll.ip arises also in the theory of propagation of surface waves 
in deep water, see [2], as an asymptotic model for small-aspect-ratio waves. 
For a; = 0, equation fll.ll) reduces to the derivative Burgers equation 

(^Uf -f- UUx')xx 0, 

whereas for b = — it reduces to the (differentiated) Vakhnenko equation 

Si] 

{ut + UUx)x + U = 0. (1.3) 

Alternatively, fll.ip with b = ^ reduces to fll.3l) after the change of variables 
{u,t) 

Equation (11.3p arises-and is known as the Vakhnenko equationin the con¬ 
text of propagation of high-frequency waves in a relaxing medium [3l [5l [6]. 
On the other hand, being written in the form 

(ti^ T Cq'^x T X (^*^) 
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it is also called the reduced Ostrovsky equation [7]: it corresponds, in the 
case /? = 0, to the equation 

CqUx CyU'llx l^'^xxx)X (^*^) 

that was derived by Ostrovsky in 1978 [8]. Therefore, it is more correct 
to name equation fll.Sp the OstrovskyCVakhnenko equation (OV), as it is 
proposed in [9]. 

Well-posedness of the Cauchy problem for the Ostrovsky equation and 
its relatives (reduced Ostrovsky equation, generalized Ostrovsky equation, 
etc) in Sobolev spaces has been widely studied in the literature, using PDE 
techniques; see [101 IHl Il2l [IS ITT] . 

On the other hand, equation fll.ip is (at least, formally) integrable: it 
possesses a Lax pair representation (see [15] appendix, (Al)) 

^xxx ^( ^xx T ^)^ (1.6a) 

i’t = ^i^xx - ui^x + Uxi^ (1.6b) 

where -0 = 0(x, t, A). 

It is because of integrable, the Ostrovsky equation can be solved by the 
inverse scattering method. Recently, Boutet de Monvel and Shepelsky for¬ 
mulate a Riemann-Hilbert problem to the initial value problem, and they 
also consider the long-time asymptotic problem, see [TB] . 

However, in many laboratory and field situations, the wave motion is ini¬ 
tiated by what corresponds to the imposition of boundary conditions rather 
than initial conditions. This naturally leads to the formulation of an initial¬ 
boundary value (IBV) problem instead of a pure initial value problem. 

In 1997, Fokas announced a new unihed approach for the analysis of 
initial-boundary value problems for linear and nonlinear integrable PDFs 
mm [T9] . The unihed method provides a generalization of the inverse 
scattering formalism from initial value to initial-boundary value problems. 
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Over the last almost two decades, it has been used to analyze boundary 
value problems for several important integrable equations with 2x2 Lax 
pairs. Recently, Lenells develop a methodology for analyzing initial-boundary 
value problems for integrable evolution equations with Lax pairs involving 
3x3 matrices [20] • Then many other integrable evolution equations with 
3x3 Lax pairs are analyzed too, see [211 SOI 123] • 

In this paper, we use Fokas and Lenells method to analyze the initial¬ 
boundary value problem for the Ostrovsky-Vakhnenko equation fll.ip on the 
half-line domain, that is, in the domain 

O = {{x,t) G M^|0 <a:<oo,0<t< T}, (1.7) 

where T < oo is a given positive constant. Assuming that a solution exists, 
we show that u{x, t) can be recovered from the initial and boundary values 
uo{x), go{t),gi{t),g 2 it) dehned by 

u{x, 0) = no(a;), 0 < x < cx), 

u{0,t) = go{t), = gi{t), u^xi0,t) = 92{t), 0<t<T. 

The main peculiarity compared with other applications of the approach of 
ca is that the Lax pair involves 3x3 matrices instead of 2 x 2 matrices. 
This difference leads to some new challenges. Apart from the 3x3 Lax 
pair, the spectral analysis of equation fll.ip on the half-line also presents 
some other peculiarities: (a) The presence of singularities in the Lax pair 
implies that it is necessary to introduce two sets of eigenfunctions. The 
eigenfunctions in the hrst set are well-behaved near z =. The eigenfunctions 
in the second set are well-behaved near z = oo. Together these two sets 
of eigenfunctions can be used to formulate a Riemann-Hilbert problem. An 
analogous situation occurs in the analysis of DP equation on the half-line 
in [21]. (b) The basic matrix eigenfunctions which are natural candidates 
for the formulation of a Riemann-Hilbert problem, are difficult to recover 
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the solution of the PDEs fl2.ip . We overcome this problem by formulating 
an associated vector Riemann-Hilbert problem, for which it is much more 
easier to recover the solution of our problem in hand, (c) The formulation 
of the Riemann-Hilber problem depends, in addition to the variables {x,t), 
on a function y{x, t) which is unknown from the point of view of the inverse 
problem. In order to obtain a Riemann-Hilbert problem whose jump matrix 
involves only known quantities, we have to reparametrize the x variable. 
This implies that we only obtain a parametric representation for the solution 
u{x, t). 

We will consider the initial-boundary value problems for fll.ip for which 
the initial and boundary values satisfy 

-Uo^^{x) + b> 0, x>0, 

-92{t) + 6>0, 0<t<T, 

as well as 

^o(^) <0, 0 <t <T. 

The assumptions in fll.9p imply that 

— Uxx{x, t) + b > 0, 0 < x < oo,0 < t < T. (1-11) 

The assumption fll.lOp is used to ensure boundedness of certain eigenfunc¬ 
tions. 

The organization of the paper is as follows. In the following section 2, we 
introduce two sets eigenfunctions. In section 3, we derive expressions for the 
jump matrices in terms of suitable spectral functions. In section 4, we derive 
residue conditions for the pole singularities of the eigenfunctions. In section 
5, we state our main result, see Theorem 15.11 


(1.9) 

( 1 . 10 ) 
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2 Spectral Analysis 


In this section, starting from the Lax pair of Ostrovsky-Vakhnenko equa¬ 
tion, see fll.hp . we define analytic eigenfunctions which are suitable for the 
formulation of a Riemann-Hilbert problem. 

Without loss of generality, in what follows we assume that 6=1. That is, 
we consider the the initial and boundary problems of the following equation 

‘^xxt Utia, “1“ SUx^xx T 0 (^*^) 


2.1 Two sets of eigenfunctions 


Let 2 : be the spectral parameter defined hj X = z^. The coefficients of the 
original Lax pair fll.bp have singularities at z = oo and also at z = 0. In 
order to have a good control on the behavior of eigenfunctions at z = oo and 
at 2 ; = 0 we introduce new forms of fll.61) . the first one appropriate at 2 ; = cxo, 
the second one at 2 ; = 0. 

Denote 


A(z) 


^ Ai(z) 0 0 ^ 

0 X2{z) 0 

0 0 A3(2;) y 


( 2 . 2 ) 


where Xj{z) = zu\ j = 1, 2, 3, here 00 


27Ti 

= 63. 
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2.1.1 Lax pair (well-controlled at z = cx)) 


Let ^ = 


^ '’P 

'Px 

y 'Pxx 


\ 




, then the Lax pair fll.hp can be written as 



^ 0 1 

0 ^ 



^x = 

0 0 

1 




z^g^ 0 

^x 

V 

u 

z-3 \ 

= 

1 

0 

—u 


^ —zuq 

1 - 

Ux 

) 


where q^{x,t) = —Uxx{x,t) + 1. 
Denote 

/ 


D{x, t) = 


\ 

( 


q 0 

0 1 

0 0 

1 1 


0 ^ 
0 

q{x,t) 

I \ 


Setting \1/ 


P{z) — Ai(z) A2 (z) A3(z) 

^ A?(z) Ai(z) ^ 

the Lax pair fl2.3p becomes 


where 


'^x - qA{z)'^ = U'^, 

+ (ngA(z) - A-i(z))§ = D'L, 


U 


^ 0 1 — 1 — 

1 — uj 0 1 — 

^^1 — 1—o; 0 j 


(2.3) 


(2.4a) 


(2.4b) 


(2.5) 


(2.6a) 
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111 

111 

111 




/ /,,2 


o;^ 0 0 


\ 


0 a; 0 
0 0 1 

(2.6b) 


/ 


Define y{x, t) by 




y{x,t) = / q{x',t'){dx' — u{x',t')dt'). 

i(o,o) 

It is well-defined, because the conservation law 


(2.7) 


qt + (uq)^ = 0. 


( 2 , 8 ) 


implies that the integral in fl2.7l) is independent of the path of integration. 
Introducing $ = we have 


l>,-[gA(z),l>] = Dl>, 
l>t + [uqA{z) - A-^(z), !>] = DI>, 

where brackets denote matrix commutator. 

We define three contours in (x, t)—domain, see Figure [H 


t* 

T 


tL 

T- 


t* 

T 






ix,t) 


O ^ 

7i 


O ^ 

72 


O ^ 

73 


(2.9) 


Figure 1: The three contours 71,72 and 73 in the (a:, f)—domain. 


And we also denote six sets which decompose the complex z—plane, see 
Figure [2] In these six sets = 1,2..., 6, the eigenvalues of A(z) and 


















Figure 2: The sets n = 1,2,..., 6, which decompose the complex 
z—plane. 


A ^{z) has the following properties. 


Di 

D2 

D3 

D, 

De 


{ReAi < ReA2 < ReAa, 
{ReAi < ReAs < ReA2, 
{ReA3 < ReAi ^ ReA2, 
{ReA3 < ReA2 < ReAi, 
{ReA2 ^ ReA3 < ReAi, 
{ReA2 < ReAi ^ RGA3, 


ReAf^ 

< 

ReA2 ^ 

< 

-^3 


ReAf^ 

< 

ReAj^ 

< 



ReAj^ 

< 

ReAh^ 

< 

-^2 


ReAg^ 

< 

ReA^^ 

< 

Ar 


ReAg ^ 

< 

ReAj^ 

< 

Ar 


ReAg ^ 

< 

ReAh^ 

< 

A3 



( 2 . 10 ) 


The solutions of fl2.9p can be constructed as solutions of the Fredholm 
integral equation 

Z e fin, 

J-y". 

‘11 

( 2 . 11 ) 
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where the contours 7 ”-, n = 1, 2, i^j = 1, 2, 3 are defined by 

{ 7 i if ReAj( 2 :) < ReAj( 2 ;) and ReA“^( 2 ;) > ReA“^( 2 ;), 

72 if ReAj(z) < ReAj( 2 ;) and ReA“^( 2 ;) < ReA“^(z), for z G r2„. 

73 if ReAi(z) > ReAj( 2 ;) 

( 2 . 12 ) 

Remark 2.1 For each n = 1,2,... ,6, the function t, z) is well-defined 
bv eauation \2.11\) for z ^ and (x.t) in the domain jg, This is because 

of the definition o/{ 7 j}i and ^2.12\) . The definition o/{ 7 j}^ implies that 

7 i : , t-t' <D, 

72 : y{x,t)-y{x',t')>D, t-t'>D, (2.13) 

73: y{x,t)-y{x',t') <11. 

The {ij)th entry of the integral equation h2.18\) involves the exponential factor 

The definition of H2.12\) implies that this factor remains bounded for z G 
when integrated along the contour 7 ”-. In fact, it is similar to that of 
^Qn{,x,t, z), except that since y{x,t) — y{x',t') can take on both signs in the 
case of'ji, the exponential is not necessarily bounded for the integration along 
7 i. However, the matrices {'y'^)ij = 'yfj for n = 1,2,.. .,6 does not involve 
integration along 71 , so we can still conclude that ^n{x,t,z) is well-defined 
in each Qn- 

Proposition 2.2 For any fixed point {x,t), ^n{x,t,z) is bounded and ana¬ 
lytic as a function of z E Qn, n = 1,2,... ,6 away from a possible discrete set 
of singularities {zj} at which the Fredholm determinant vanishes. Moreover, 
^n{x,t,z) admits a bounded and continuous extension to Qn o^nd 

^n{x,t, z) = I0{-), Z —>■ OC, Z e Qn- (2-14) 

"Z 
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2.1.2 Second Lax pair (well-controlled at z = 0) 

Setting \&o = then the Lax pair fl2.3p becomes 

'I'Oo; - A( 2 :)'Lo = 

= Vo^o, 

where 


(2.15) 


Un = 


ZUr, 


u 

0 

0 ’ 


! 1 

1 

1 

0 

a;2 

0 


1 

1 

1 

1 0 

0 

1 > 


1 1 

1 

1 > 


(2.16a) 


/ 


ILcr 


Va = -^ 


0 


1 — U 1 — U) 


\ 


\ 


1 — u 0 1 — o;^ 

1 — 1 — a; 0 


-zu 


V 


a; 0 0 

0 ^2 0 

0 0 1 


u^x 



(2.17) 


Setting <l)o = then the Lax pair of <ho is 

l>ox- [A(;2),<ho] = 7/o$o, 
l>ot-[A-H;^),<ho]=H$o, 

whose solutions can be constructed as solutions of the Fredholm integral 
equation 

(l>on(x,t,z))p = 5 ,,+ [ ^ G ff. 

' ‘i'J 

( 2 . 18 ) 

where q”- are defined as fl 2 . 12 p . 

Remark 2.3 For each n = 1,2,... ,6, the function <Fon(a^, t, z) is well-defined 
bv eauation h2.13i) for z ^ and (x.t) in the domain jri' This is because 
of the definition o/{7j}f and \2.12\) . The definition of {jj}l implies that 


7 i : a; — x' > 0, f — F < 0, 

72 : X — x' >D,t — t'>D, (2-19) 

73 : a; — x' < 0 . 
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The {ij)th entry of the integral equation H2.18\) involves the exponential factor 


^{Xi(z)-Xj{z))(x-x')+{X^ ^(z)-Xj l(z))(t-t') 


The definition of ^2.12\) implies that this factor remains bounded for z G fin 
when integrated along the contour 7 ”-. 


Proposition 2.4 For any fixed point {x,t), ^on{,x,t,z) is bounded and ana¬ 
lytic as a function of z E fin, n = 1,2,... ,6 away from a possible discrete set 
of singularities {zj} at which the Fredholm determinant vanishes. Moreover, 
^on{x,t,z) admits a bounded and continuous extension to fin and 


t, z) — I + ^$0 ^ ^ “1“ 0{z^), Z —>■ 0, Z E. fin, (2.20) 


where 


i" = -5«J, if = 


( 2 . 21 ) 




with r = 


OJ Ul OJ 


\ 


b? LkP' 


111 


and f = [A, r], A = diag{u, 1}. 


2.1.3 Further properties of ^n and <l)on 

Now, noticing that <i>„ and <l>on are related to the same linear system of 
PDFs fl2.3p . tracing back the way that the differential equations for and 
•hon were derived, leads to the following proposition. 


Proposition 2.5 The functions <!>„ and <hon are related as follows, 

^nix,t,z) = p-\z)D-\x,t)Piz)%n{x,t,z)e^^-y^^'^^^^^^\ ( 2 . 22 ) 


where P{z) and D{x,t) defined as ^2.4\), and y{x,t) is defined as ^2. 7|). 
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Noticing that it is a straight calculus to show P~^{z)D~^{x,t)P{z) is inde¬ 
pendent of 2 ;. 

From fl2.14p . fl2.20p and fl2.22p we derive the following expansion of <Fn(x, f, z) 
as z 0, 


^n(x,t, z) 


p-\z)D-\x, t)P(z) (^I -t- z |-f F + (x - 2 /)a| 
+^2 |_|f - f(x - y)Tt + + 0(z3)j 


(2.23) 


2.2 Symmetries 

Proposition 2.6 z) satisfies the symmetry relations: 


$„(x, t, ujz) = A^n{x, t, z)A ^, with A = 


/ 0 1 0 \ 
0 0 1 
yl 0 Oy 


<Fn(x, t, z) = B^nix, t, z)B with B = 


/ 0 1 0 \ 
1 0 0 
yo 0 ly 


= B 


-1 


^nix,t, z) = C^n{x,t,uj'^z)C ^, with C = 


$„(x, t, z) = V^n(x, t, uz)V with V = 


/ 0 0 1 \ 
0 1 0 
1 0 0 

/ 1 0 0 \ 
0 0 1 
yo 1 Oy 


= C 


-1 


= v 


-1 


The functions <Fon(x,f, z) satisfies the same symmetry relations. 
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3 Sectionally Meromorphic Functions 


Let {‘Ln}i and {$on}i denote the eigenfunctions defined in Section 2. We 


have well control over near 2 : = oo. On the other hand, we have well 
control over <hon near z = 0. Therefore, we will introduce a radius R > 0 and 
formulate a Riemann-Hilbert problem by using the for \z\ > R and the 
$on for \z\ < R. 

Dehne sets {Dn}\^ by (see Figure [3]) 



Figure 3: The sets n = 1,2, which decompose the complex 

2 ;—plane. 


= On n {| 2 :| > R}, n = 1,2 ,... ,6, 
Dn+e = r) {\z\ < R}, n = l,2, 


(3.1) 


Since the map F : (x, t) {y,t),y = y{x, t) is a bijection from the domain 
fll.71) onto F{Q), we can dehne functions {Mn{y,t,z)}\^ for {y,t) G F{Q) by 



z e Dn,n = 1,2,... ,Q, 
z E Dn, n = 7, 8 ,..., 12. 


(3.2) 
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The Mn defined in fl3.2p are bounded and analytic on the complex 2 ;—plane. 
Reminding that the relation condition between the two eigenfunctions 
and <hon, we need formulating the Riemann-Hilbert problem in terms of the 
row vectors z/„ defined by 

1 1 j z), z e Dn, n = 1,2,.. .12. (3.3) 

Let M and v denote the sectionally meromorphic functions on the com¬ 
plex z—plane which equal and Vn respectively for z G Dn- 


Lemma 3.1 The function M obeys the symmetries 

/ 0 1 0 \ 

• M(x, f, ojz) = AM{x, t, z)A~^, with A = 


V 


0 0 1 
1 0 0 


/ 


M{x, f, z) = BM{x, t, z)B , with B = 


0 1 0 
1 0 0 
0 0 1 




= B-\ 


( 


M{x,t^ z) = CM{x,t,(z‘^z)C ^, with C = 


M{x, t, z) = 'DM{x, t, uz)V \ with V = 


V 


0 0 1 
0 1 0 
1 0 0 

1 0 0 
0 0 1 
0 1 0 




= C-h 




/ 


= V 


-1 


Proof: This is a consequence of equation fl3.2l) and the symmetry properties 
of the <1)„ and the $on- Noticing that 


P(z) = P{ujz){A), P{z) = P{z)B. (3.4) 

And the other two symmetries can be obtained by these two symmetries. □ 
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3.1 The jump matrices 

We define spectral functions Sn{z) by 

Sn{z) = Mni0,0,z), zeDn, n = (3.5) 

The tracelessness of the matrices {17, V} and {Uq, Vq} implies that 

det5„(;2) = l, n = 1,2, (3.6) 

The exponential factor on the right-hand side of fl3.2p has been 

included because it ensures that the jump matrices introduced in the next 
proposition depend on x only through the function y{x,t). 

Proposition 3.2 For each n = 1, 2,..., 12, the function Vn is bounded and 
analytic in (away from the possibly empty discrete set {zj}). Moreover, 
each Un has a continuous and bounded extension to Dn- The function v 
satisfies the jump conditions 

TZn = T^mJm,n, Z E Dm, U, m = 1, . . . , 12, H ^ m, (3.7) 

where the jump matrix Jm,n(y, t, z), Jm,n = i^ defined by 

Jm,n = n,m G (1 ,..., 12}. (3.8) 

Proof: The analyticity and boundedness properties of the Un follow from 
the properties of the and the <hon established in Section 2. 

From the dehnition of (12.111) and $on (I2.18p . we can deduce that 

Mn = (3.9) 

where J(z) is a matrix independent of (x,t). 

Evaluation at a; = f = 0 yields J = SjjfSn- Multiplying fl3.9p by 
^ 1 1 1 j from the left, we obtain the jump condition fl3.7p with 
given by □ 


16 












4 Residue Conditions 


If the Un have pole singularities at some points {zj},Zj G C, the Riemann- 
Hilbert problem needs to include the residue conditions at these points. We 
will assume that all Zj lie in the interiors of the sets {Dn}f, singularities in 
the interiors of the sets {Dn\Y t)e avoided by choosing R large enough. 
The residue conditions can be found by relating the to another set of 
solutions of fl2.15p . denoted by {^j}i , which are dehned by 


Hj{x,t,z) 


1+ '^^\Udx'+Vdt'){x',t')fij{x',t',z), 


'7j 


where {7j}f are contours shown in Figured! 


J = 1,2,3, 
(4.1) 


4.1 A matrix factorization problem 

Let us dehne the 3 x 3—matrix value spectral functions 5 ( 2 ;) and S{z) by 



^ 3 {x,t,z) = 

(4.2a) 


fj.i{x,t,z) = fj, 2 {x,t, 

(4.2b) 

Thus, 

s(^) = h 3 ( 0 , 0 , 2 ), S'( 2 :) = ^ 1 ( 0 , 0 , 2 :). 

(4.3) 

Lemma 4.1 Due to the symmetries of Mn, see Lemma\d.l\ we just need to 
calculate Si. The Si defined in (|,?. 5\) can he expressed in terms of the entries 
of s{z) and S{z) as follows: 


/ Sii 0 0 \ 

Si- S21 0 

m23is) 1 

\ 31 m33(s) / 

(4.4) 

where 

mij denote that the {i,j)—th minor of s. 
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Proof: Let 7 ^° denote the contour (XqiO) {x,t) in the (x, t)— plane, here 
Xq > 0 is a constant. We introduce z; Xq) as the solution of fl4.ip 

with j = 3 and with the contour 73 replaced by 7 ^°. Similarly, we dehne 
Mi{x,t, z] Xo) as the solution of fl3.2p with 73 replaced by 7 ^°. We will 
hrst derive expression for 5'i(z;Xo) = Mi(0,0,z;Wo) in terms of S{z) and 
5 ( 2 ;; Wq) = /i3(0,0, z; Xq). Then fl4.4p will follow by taking the limit Xq —)■ 00 . 

First, from the dehnition of Mn fl3.2p . we have the following relations: 

{ Mi{y, t, z; Xo) = //i(x, t, Ri{z; 

Mi{y, t, z] Xq) = /i 2 (x, t, 

Mi{y,t,z]Xo) = 

(4.5) 

Then we get i?i( 2 ;;Xo) and Ti(z;Xo) are dehned as follows: 

R,{z;Xo) = (4.6a) 

T,{z-Xo) = e-^(^)^°M„(2/(Xo,0),0,^;Xo)e"(^°’°)^(^\ (4.6b) 

The relations fl4.5p imply that 

s(^;Xo) = S,{z-Xo)T,-\z-Xo), S{z) = S,{z-Xo)R^\z-Xo). (4.7) 

These equations constitute a matrix factorization problem which, given { 5 ( 2 ;), 5'(2:)} 
can be solved for the {i?i. S'!, Ti}. Indeed, the integral equations fl3.2p to¬ 
gether with the dehnitions of {i?i, 5'i,Ti} imply that 

f {RiU-,Xa))ij = 0 if 7,^ = 7i, 

I (Si(3;X„))., =0 if 7 ^ = 72, (4.8) 

( (T,(2iV„))., =0 if 7.(j = 73. 

It follows that fl4.7p are 18 scalar equations for 18 unknowns. By computing 
the explicit solution of this algebraic system, we hnd that {S'i(z;Xo)} are 
given by the equation obtained from fl4.4p by replacing {S'i( 2 ;), s(; 2 )} with 
{S'i(z; Xq), s(z; Xq)}. Taking Xq —)■ cx) in this equation, we arrive at fl4.4p . 

□ 
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4.2 The residue conditions 


Since /i 2 is an entire function, it follows that M can only have singularities 
at the points where the have singularities. We infer from the explicit 
formulas fl4.4p that the possible singularities of M are as follows: 

• [M ]2 could have poles in Di at the zeros of sii( 2 ;); 

• [Mjs could have poles in Di at the zeros of 

We denote the above possible zeros by {zj}^ and assume they satisfy the 
following assumption. 

Assumption 4.2 We assume that 

• sii( 2 ;) has ni possible simple zeros in Di denoted by {zj}i^; 

• m^.z^z) has N — ni possible simple zeros in Di denoted by 

and that none of these zeros coincide. Moreover, we assume that none of 
these functions have zeros on the boundaries of the Dn’s. 

We determine the residue conditions at these zeros in the following: 

Proposition 4.3 Let {Mi} be the eigenfunctions defined by H3.W and as¬ 
sume that the set {zj}i of singularities are as the above assumption. Then 
the following residue conditions hold: 



(4.9a) 


1 < j < ni, Zj G Di 



e023bd[M(zj)]2, 


(4.9b) 


< i < N, Zj G Di 


where f = ^, and 9ij is defined by 

%(z) = (A, - Xj)y+{X-^ - X-% t,j = 1,2,3. (4.10) 
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Proof: We will prove fl4.9ap .From the relation 




(4.11) 


For i, j = 1, 2, 3, let Oij = (A, — Aj)a;+ (Aj ^ — A^ ^)t. In view of the expressions 
for Fi given in fl4.4l) . the three columns of fl4.1ip read: 


[Ml] 26^1 
[Mijge^i^^-^) 


[fJ‘2]lSll{z) + [fj.2]2e^^^ S2l{z) + [fi 2 ] 3 e^'^^ S 3 i{z) , 


+ Nse»” 




Sll 

1 


2 3 


mssis)' 


^ 23 ( 5 ) 

Sll 


(4.12a) 

(4.12b) 

(4.12c) 


In order to prove fl4.9ap . we suppose that Zj G Hi is a simple zero of 
511 ( 2 ;). Solving (I4.12ap and (I4.12cp for [^ 2)2 and [/i 2]3 and substituting the 
result in to fl4.12bp . we hnd 



Sii 521 "S21 S21 


Taking the residue of this equation at Zj, we hnd the condition fl4.9ap in the 
case when Zj G Hi. Similarly, we can get the equation fl4.9bp . □ 


5 The Riemann-Hilbert problem 

The sectionally analytic function p{y, t, z) dehned in section 3 satishes a 
Riemann-Hilbert problem which can be formulated in terms of the initial 
and boundary values of u{x,t). By solving this Riemann-Hilbert problem, 
the solution of fl2.ip can be recovered in parametric form. 

Theorem 5.1 Suppose that u{x,t) are a solution of ^2.1\) on the half-line 
domain H ( [i. 7| ) with sufficient smoothness and decays as x ^ 00 . Suppose 
that the initial and boundary values {uo{x), go{t), gi{t), g 2 it)} defined in (11) 
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satisfy the assumptions jn). Then u{x, t) can be reconstructed from the 
initial value and boundary values as follows. 

Use the initial and boundary data to define {$n(0,0, z)}® and {$on(0,0, z)}i 
via the integral equations h2.11\) and H2.18\) . respectively. Define spectral func¬ 
tions Sn{z), n = 1,... ,12, by 

Sn{z) = %n-6{0,0,z), Z E D^, Tl = 1, . . . , 6, 

S'„( 2 ;) = P"^(2;)D(0,0)P(2;)I>„(0,0,2;), z E D^, n = 7,..., 12, 

where P{z) and D{x,t) are defined in Iji2.4^ . Define the jump matrices 
Jm,n{y,t,z) in terms of the Sn by equation H3.8\) . Define the spectral s{z) 
and S{z) by equation \4.2^ . Assume that the possible zeros {zj}i of the 
functions Sn( 2 ;), 77133 ( 5 ( 2 ;)) are as in assumption \4.2\ 

Then the solution u{x, t) is given in parametric form by 

u{x,t) = u{y{x,t),t), (5.2) 

where 

t) = y + \imj- {y, t, ^) - l) , (5.3a) 

«(,.*) = (5.3b) 

with i'A){ii^t,z) is the third column of row-vector value function i'{y,t,z) 
which satisfies the following 3x3 vector Riemann-Hilbert problem: 


• v{y, t, z) is sectionally meromorphic on the complex z—plane with jumps 
across the contour D^ n n,m = 1, 2,..., 12, see Fiqure\^ 


• Across the contour Dn H Dm, 


V satisfies the jump condition \3. 7|) 


• v{y,t,z) = 1 lj+0(i), z^-cx). 

• has simple poles at z = Zj, 1 < j < ni. has simple poles at 
z = Zj, ni + 1 < j < N. The associated residue condition is showed in 
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the following: 


ReSz=z = ’^33(^(2j)) p9i2{zj)^{i)(^^ ^ \ 

^ ^ Sll(zj)s2i{zi) ^ 

^ < 3 <ni,Zj e Di 

ReSz=zi^^'^\yitiZ) = — «ii(zj) 

^ ^ ^ niz3{s{zj))m2i{s{zj)) ^ ^ (5.4b) 

ni < 3 < N, Zj G Di 


where f = ^, and 9ij is defined by 


9ij{z) = {Xi - Xj)y + {X- - Xj )t, i,j = l,2,3. (5.5) 


For each zero Zj G Di, there are five additional points, 


uzj, oj'^Zj, kj, ujZj, u^Zj, 


at which v also has simple poles. The associated residues satisfy the 
residue conditions obtained from and the symmetries of Lemma 

mi 

Proof: The residue conditions (15.dh are obtained by multiplying the con¬ 
ditions in fl4.9p by ^ 1 1 1 ^ from the left. It remains to prove fl5.3p . 
Equation fl2.23p implies that 

u(y,t,z) = (^ ui 102 U3^+0{z^), as z^O, (5.6) 

where 

ui= (^1 + zuj{x -y) + ^z^u‘^{x - yf^ , (5.7a) 

V 2 = (^1 + zoj‘^{x -y) + ^z^uj{x - yf^ , (5.7b) 

V 3 = (^ +z{x-y) + ^z^{x-yf^ , (5.7c) 

Thus, from the equation fl5.7cp we can get fl5.3p . □ 
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